Intrinsic and extrinsic corrugation of monolayer graphene deposited on Si0 2 
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Using scanning tunneling microscopy (STM) in ultra high vacuum and atomic force microscopy, 
we investigate the corrugation of graphene flakes deposited by exfoliation on a Si/SiCb (300 nm) 
surface. While the corrugation on S1O2 is long-range with a correlation length of about 25 nm, 
some of the graphene monolayers exhibit an additional corrugation with a preferential wave length 
of about 15 nm. A detailed analysis shows that the long range corrugation of the substrate is also 
visible on graphene, but with a reduced amplitude, leading to the conclusion that the graphene 
is partly freely suspended between hills of the substrate. Thus, the intrinsic rippling observed 
previously on artificially suspended graphene can exist as well, if graphene is deposited on SiCb- 
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Since it was believed, based on the Mermin- Wagner 
theorem, that two-dimensional (2D) crystals are not sta- 
ble at finite temperature [1 j, it came as a surprise that 
monolayer graphene could be stabilized on a Si/Si02 
substrate [2j 02 [4]. Due to its peculiar properties like, 
e.g., a linear dispersion leading to Klein tunneling [5], 
high room-temperature mobility allowing quantum Hall 
steps at 300 K [6], a low spin-orbit interaction benefi- 
cial for spintronic devices [7] or tunable spin-polarized 
edge states [8], graphene studies have become a ma- 
jor issue in solid-state physics. Already the first trans- 
port results not exhibiting weak localization lead to the 
speculation of a curved surface [9] acting as a phase- 
breaking field [10 j. Such a curvature was indeed observed 
by microscopic electron diffraction of suspended mono- 
layer graphene sheets [llj. The lateral wavelength of the 
isotropic curvature is estimated to be A = 10—25 nm with 
an amplitude of A ~ 1 nm. The rippling has been re- 
produced theoretically by Monte-Carlo simulations with 
a preferential wavelength of 8 nm barely depending on 
temperature [12] . It is argued that the anharmonic cou- 
pling between bending and stretching modes in graphene 
causes the rippling and is responsible for the stability of 
the 2D crystal. Moreover, it has been shown experimen- 
tally that rippling can ultimately limit the mobility of 
graphene at 300 K, if defects are avoided [13 

However, the corrugation properties of graphene de- 
posited on a substrate as typically used in transport ex- 
periments are not clarified. In particular, previous scan- 
ning probe studies [14] [15] [16] [17] revealed only corruga- 
tions, which are attributed to the roughness of the under- 
lying substrate, except for a torn flake manipulated by 
AFM [16 . Thus, it is unclear, if the intrinsic tendency for 
rippling persists on the substrate, thereby ultimately lim- 
iting mobility and influencing weak localization, which 
has meanwhile been found, at least, for some samples 
[16]. Here, we demonstrate that a regular short wave- 



length corrugation with A ^ 15 nm and A c± 1 nm, which 
is not induced by the substrate, can even prevail the sub- 
strate corrugation. Since A and A are close to the values 
found on suspended graphene [11], we attribute this ad- 
ditional corrugation to the intrinsic rippling of graphene 
[T2] . Moreover, we find that the long-range corrugation 
implied by the substrate has lower amplitude on graphene 
than on S1O2 suggesting that our high-mobility graphene 
is partly suspended between hills of the substrate, which 
might favor the development of the intrinsic short-scale 
rippling. 

The graphene sample is fabricated by the mechanical 
exfoliation technique as described in [2j [18]. Using an 
optical microscope, a graphene flake containing a mono- 
layer region is identified. Raman spectroscopy is used 
to confirm the number of layers with the help of the 2D 
line [19]. Fig. [TJa) shows three spectra measured in dif- 
ferent areas of the flake visible in Fig. TJb). In addition 
to a large monolayer area, we find smaller bi- and multi- 
layer areas as marked by 1L, 2L and MuL, respectively. 
Next, gold contacts with a 10 nm Cr seed layer were de- 
posited and structured with a lift-off process. Except for 
one side, the graphene sample is completely surrounded 
by the gold electrode as shown in Fig. [TJb). The mo- 
bility of graphene samples prepared identically but with 
several contacts is \i — 1 — 1.5 m 2 /(Vs) at 300 K [20] . 
i.e. comparable with high-mobility values obtained by 
other groups [13] [21]. In order to remove residual re- 
sist and adsorbates, the sample was rinsed in isopropanol 
and acetone, baked out to 150° C in air for four hours 
and, additionally, baked out to the same temperature 
in ultrahigh-vacuum (UHV) for three hours. The gold 
film served as the electrical contact for scanning tunnel- 
ing microscopy (STM) and was used to prepare the STM 
tip (etched W) by voltage pulses. In order to find the 
graphene within the UHV- STM, we used an optical mi- 
croscope with a focal length of 30 cm and 5 /im lateral 
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FIG. 1: (Color online) (a) Raman spectra using laser light of 
wavelength 532.1 nm. The spectra are recorded on different 
areas of the graphene flake and the attribution to monolayer, 
bilayer and multilayer is marked, (b) Optical image of the 
graphene sample with gold contacts. The different areas of 
monolayer (1L), bilayer (2L) and multilayer (MuL) as iden- 
tified by Raman spectroscopy are indicated. Insets show two 
STM images acquired at the edge of the gold contact, which 
are used for orientation, (c) Microscopic image of the STM 
tip in tunneling contact with the graphene in UHV recorded 
with a long-range optical microscope. 

resolution. Fig.[TJc) shows a microscopic image of the tip 
approaching the graphene flake. The home-built STM 
features an xy stage for lateral positioning and operates 
in UHV (2-l(T 7 Pa) as described elsewhere [22j|23]. The 
xy stage has been checked to move accurately within 10 
%. Topographic images are recorded applying a bias V 
to the tip. Spectroscopic dl/dV curves are measured by 
lock-in technique using a modulation voltage V mo d af- 
ter stabilizing the tip at voltage V^tab and current i s t a b- 
STM images obtained at the intentionally irregular edge 
of the gold contact are displayed in the insets of Fig.[TJb). 
They are used for orientation by STM and aid in finding 
monolayer, bilayer and multilayer areas on the flake. 

Atomic force microscopy images are taken under ambi- 
ent conditions in tapping mode using either a commercial 
cantilever with a silicon tip or an ultrasharp tungsten tip 
with a tip radius of 1 nm, which is attached at the bot- 
tom of a Si cantilever [24] . Imaging has been carried out 
in the attractive regime [25] with an oscillation frequency 
slightly above resonance. 

Fig. [2^b) and (c) show two atomically resolved STM 
images recorded on the monolayer (1L) and the multi- 
layer (MuL) of graphene, respectively. One observes a 
hexagonal pattern on the monolayer and a triangular pat- 
tern on the multilayer in accordance with previous results 
P~4j[T7]. The atomic resolution is still perceptible at lower 
resolution in Fig.[2|a) and (d), but an additional irregular 
corrugation appears. Line sections shown in Fig. [2|e) re- 
veal that the corrugation height on this monolayer region 



is about 0.6 nm, a factor of three larger than on the mul- 
tilayer region. Notice that corrugation heights pretended 
by charged defects are typically smaller by more than an 
order of magnitude [26]. A representative dl / dV / {I /V) 
curve of the monolayer, which is known to represent the 
local density of states (LDOS) [27 , is shown in Fig. [2^f ) . 
As expected for graphene, it shows a V-like shape with 
the Dirac point at about V = —20 mV. It does not change 
significantly across the area depicted in Fig. J2jd). Note 
that graphene on SiC(OOOl) exhibits a more complicated 
dl/dV spectrum [28j [29j [30] . Note also, that we did not 
observe the phonon gap found at 4 K recently [17], al- 
though we used several different micro- and macrotips. 
This is tentatively attributed to the different tempera- 
ture of the two experiments. The good representation of 
the graphene LDOS by dI/dV/(I/V) and the continu- 
ous atomic resolution across Fig. |5|d) demonstrates the 
absence of resist on this part of the surface. Indeed, we 
found several large areas without resist, but still also ar- 
eas where a remaining resist is apparent within the STM 
images. The fact that the Dirac point is observed close to 
V shows a negligible influence of charging adsorbates, 
which are probably removed in UHV. 

Fig.[3^a) shows a large area STM image of the graphene 
monolayer. One observes a rather regular corrugation 
with amplitudes of more than 1 nm and a preferential dis- 
tance between hills of about 15 nm. The histogram [32] is 
Gaussian with a full width at half maximum (FWHM) of 
0.78 nm. The rms roughness is determined to be 0.36 nm. 
The latter two values fluctuate across the monolayer area 
by about 25 % (average rms roughness 0.32 nm), but the 
corrugation length scale remains constant. For compar- 
ison, Fig. [3jc) shows an AFM image of the bare Si02 
surface recorded in tapping-mode. It is obvious that 
the corrugation on the substrate is lower and exhibits 
a larger length scale than the corrugation on graphene. 
The FWHM of the histogram of Fig. [3jc) is 0.48 nm 
and the rms roughness is 0.22 nm (average of all data 
(0.25 ± 0.05) nm). We used four different tips, including 
two ultrasharp tips with a tip radius of 1 nm [24 , but 
always observed the same corrugation height and length 
scale, which, in addition, is in good agreement with pre- 
vious results [T5] . 

Careful inspection of the graphene monolayer in 
Fig.^a) and (b) shows a modulation on top of the short- 
range corrugation exhibiting a similar length scale as the 
corrugation on the Si02 substrate. In order to disen- 
tangle the two contributions, we used two-dimensional 
autocorrelation functions of the images (a) and (c) [31], 
which are shown in Fig. [4ja) and (b). Note that both 
figures use the same color scale, i. e. the intensities are 
directly comparable. A long-range structure represented 
by a central spot surrounded by four bright areas is vis- 
ible in both images although with weaker intensity for 
graphene (a). Fig. |4jc) shows the autocorrelation func- 
tion of the graphene after removing the long-range part 
by high-pass filtering using a smooth (first order Butter- 
worth) cutoff at wavelength A = 20 nm. A preferential 
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FIG. 2: (Color online) (a) Constant current STM image of multilayer graphene (0.4 V, 258 pA) (raw data), (b) Higher 
resolution STM image of the same area as in (a) (0.4 V, 198 pA). (c) High resolution STM image of graphene monolayer area 
(IV, 394 pA). (d) Lower resolution STM image of monolayer graphene (0.5 V, 292 pA) (raw data), (e) Line section along the 
lines depicted in (a)(MuL) and (d)(lL). (f) Normalized dl/dV spectrum of monolayer graphene: V^tab^O.TV, 7 s tab=300pA, 
Knod=20mV. 
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FIG. 3: (Color online) (a), (b) 3D and 2D constant current STM image of monolayer graphene (IV, 207 pA). (c) 3D tapping- 
mode AFM image of the S1O2 substrate (resonance frequency 326.4 kHz, force constant 47 N/m, excitation frequency 326.5 kHz, 
oscillation amplitude 18 nm, constant amplitude feedback, setpoint 90%). 



short-range distance is clearly visible by the strong spots 
around the center, but there is additional multiple corre- 
lation albeit not with high symmetry. Notice that creep 
and drift effects are carefully removed from the STM- 
image using the atomic resolution. We checked the pref- 
erential directions visible in both correlation patterns. 
They are rotating arbitrarily across the Si02, but we 
find indication for a preferential orientation on graphene 
[32] . Fig. [4^d) shows radial line sections averaged over 
all angles of the correlation functions. The Si02 (black 
curve) exhibits a correlation length of about 25 nm com- 
parable to [15] and, in addition, a very weak maximum 
at about 50 nm indicating a slight preferential distance 
between hills. The gray curve shows the radial line sec- 
tion of the high-pass filtered image of graphene visible in 
Fig. |4jc). It exhibits a damped oscillation with a wave 
length of about 14 nm looking very similar to correla- 
tion functions of liquids. For comparison, the radial line 



section of the autocorrelation function of the AFM im- 
age after applying the same high-pass filter is shown as 
a dotted line and does not exhibit any structure. We 
observe the damped oscillation also after high-pass filter- 
ing the original STM image in Fig. [3]^a) and determining 
the autocorrelation function afterwards. In addition, we 
checked that the obtained wave length is quite robust 
with respect to slight parameter changes of the filtering 
process. Depending on the filtering procedure and the 
selection of the minima or maxima used for determining 
the wave length, A is varying by at maximum ±2 nm 
with a mean value of A = 15 nm. 

From these results, we conclude that monolayer 
graphene can exhibit a rather regular, liquid-like short- 
scale corrugation not induced by the substrate. It is 
similar in height and wavelength to the one observed on 
suspended graphene [11 and calculated for freely sus- 
pended graphene by atomic Monte-Carlo methods [T2] . 
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FIG. 4: (Color online) (a) 2D autocorrelation function of the 
drift corrected STM image of monolayer graphene as shown in 
Fig. [3^ a), (b) 2D autocorrelation function of the AFM image 
of the Si02 surface as shown in Fig.[3|c), same color scale as in 
(a), (c) FFT high pass filtered image (1. order Butterworth, 
20 nm) of the 2D autocorrelation function shown in (a), (d) 
Radial line sections of the unfiltered AFM autocorrelation im- 
age (b) (black) , of the high pass filtered STM autocorrelation 
image (c) (gray), and of the high pass filtered AFM autocor- 
relation image (dotted). The inset shows a larger view of the 
area marked by the dashed line. The black and gray graphs 
are offset vertically for clarity. 



This observation suggests that the graphene might be 
partly suspended even on the Si02- In order to verify 
this, we compared the Fourier components of the images 
of graphene and Si02 within the large wavelengths re- 
gion A = 40 — 60 nm [32] . We indeed find that graphene 
exhibits an up to 40 % lower amplitude at all wave- 
lengths within this wavelength range. This implies that 
the graphene does not follow the substrate corrugations 
exactly, but instead includes areas not in contact with 
the substrate. We believe that this partly free-standing 
configuration is the prerequisite for the intrinsic rippling 
of graphene. Indeed, we have also found other graphene 
flakes prepared nominally identically by the same person 
which did not show the intrinsic rippling. This high- 
lights that tiny details of the preparation process can 
have large impacts on the physical properties, maybe ex- 
plaining the contradicting weak localization properties in 
different studies [9j [16] . A comparison of Raman data of 
several samples indicates that samples with intrinsic rip- 
pling are more frequent than samples without [32 . Fu- 
ture systematic studies using our method to determine 
the intrinsic rippling are required. 

In summary, we have investigated high-mobility 
graphene flakes on Si02 by UHV-STM and scanning tun- 
neling spectroscopy. A mesoscopic corrugation decreas- 
ing in amplitude with increasing thickness not induced by 
the substrate is identified. It is rather regular exhibiting 
liquid-like correlation properties with a preferential wave- 
length of about 15 nm and a rms roughness of 0.32 nm. 
Since the long-range corrugation of the substrate is ad- 
ditionally visible on graphene, but with smaller ampli- 
tude than on the substrate, we conclude that the rippled 



FIG. 5: (Color online) Height histograms of SiCb (blue) and 
graphene (red) corresponding to the images in Fig. 3 of the 
main text. On average, the graphene histogram exhibits a 
larger FWHM value compared to SiCb. 

graphene is partly freely suspended. 

We gratefully acknowledge useful discussions with G. 
Giintherodt, U. D. Schwarz and H. Kurz and financial 
support by FOR 912, project TP 6 of the Deutsche 
Forschungsgemeinschaft and by the German Federal Min- 
istry of Education and Research (BMBF) under contract 
number NKNF 03X5508 ("ALEGRA"). 

Supplementary information 

Histogram of Corrugation 

Figure [5] shows representative histograms of the height 
values on the Si02 substrate and the graphene sample as 
obtained by atomic force microscopy (AFM) and scan- 
ning tunneling microscopy (STM), respectively. The his- 
tograms are based on images of 200 x 200 nm 2 including 
512 x 512 measurement points each. These images are 
shown in Fig. 3 of the main text. The histogram curves 
are normalized in order to enclose the same integral area. 
Gaussian fits are added as guides to the eye. Obviously, 
the height histogram of Graphene is broader exhibiting a 
full width at half maximum (FWHM) of 0.78 nm, while 
Si0 2 exhibits a FWHM of only 0.48 nm. While the latter 
value is quite robust across the surface, the former value 
is changing by about 25 % across the surface, but is on 
average larger than the FWHM on Si02- This demon- 
strates that the rippling of graphene is not limited to 
substrate corrugations. 



Tip Dependence of Substrate Corrugation 

Figure [6] shows two AFM images of the Si02 substrate 
obtained in tapping mode by two different tips. Rep- 
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resent ative scanning electron microscopy (SEM) images 
of the tips provided by the supplier are added as in- 
sets. They are cross-checked by SEM images after the 
measurements. These images are in agreement with the 
supplier images, but due to a moderate background pres- 
sure, the ultrasharp W-tip fastly burnt away during SEM 
imaging. The left AFM image is measured by a cantilever 
with a Si tip having an apex radius of about 10 nm, while 
the right AFM image is measured with a W tip of nomi- 
nal tip apex radius of 1 nm mounted at the bottom of the 
Si tip. Both images exhibit a very similar lateral scale 
of the corrugation. The rms values are 220 pm and 270 
pm, respectively, showing that the tip properties do not 
significantly change the observed corrugation of the sub- 
strate. We cross-checked this result using more than ten 
images obtained with both tips, which lead to fluctua- 
tions in rms values of less than ± 10%. 



Comparison of Correlation Functions 

Figure [7] shows a comparison of the correlation func- 
tions of the topographic images of graphene (top row) 
and Si02 (bottom row). The left column shows the cor- 
relation functions directly obtained from the topography. 
The middle and right column show the high-frequency 
and low-frequency part of the correlation function sep- 
arated by a smooth first-order butterworth filter with 
cutoff wavelength A = 20 nm. While the large wave- 
lengths visible in the right columns are quite similar in 
both images, the rather regular short wavelength pattern 
(middle column) is only observed on graphene. Note that 
correlation functions within the same columns are scaled 
to the same contrast, i. e., their brightness is directly 
comparable. 



Reproducibility of Correlation Functions 

Figure [8] shows three images (top row) obtained in dif- 
ferent areas of the monolayer graphene, separated by sev- 
eral /im. They exhibit different rms values of 270 pm, 220 
pm and 360 pm (left to right), respectively. The mid- 
dle row shows the corresponding correlation functions. 
The bottom row shows the same correlation functions 
after high-pass filtering using the identical procedure as 
in Figure 3. Notice that the scales of the images are dif- 
ferent. While the long-range corrugation (middle row) is 
oriented in different directions, the short range corruga- 
tions (bottom row) exhibit the same wave length of about 
15 nm in all three areas and show a similar orientation 
as marked by the white lines. The latter indicates that 
the the short-scale corrugation is related to the atomic 
directions in graphene. But since the directions of the 
rippling are fluctuating slightly across the sample sim- 
ilar to the orientations of molecules in liquid crystals, 
more data are required in order to substantiate the rela- 
tion between atomic directions and preferred directions 



of rippling. 

Figure [9] shows topographic images of different areas 
of the Si02 substrate (top row) and the corresponding 
correlation functions (bottom row) . They exhibit similar 
length scales of corrugation but clearly different orienta- 
tions, which might be due to the polishing procedure of 
the sample. 



Wave Length Dependence of Corrugations 

Figure [To] shows the radially averaged Fourrier trans- 
forms (FT) of the topographic images of graphene (black 
curve) and Si02 (grey curve). The FT curves are based 
on images of the same lateral size and resolution. The 
curves are averaged using 10 images of 200 x 200 nm 2 
from different areas of the surfaces. Obviously, the corru- 
gation of Si02 is larger at large wave length (small wave 
number) down to 30 nm by up to 66 %, while the corru- 
gation on the graphene is larger at smaller wavelengths 
(larger wave numbers) by up to a factor of 2.6. The 
individual FTs of different Si02 areas are nearly indis- 
tinguishable, i.e they are not deviating from the width of 
the grey curve in Fig. 5. However, the FTs of different 
graphene areas are shifted vertically due to the fluctu- 
ating rms value of corrugation, which leads to crossing 
points of the two curves varying between wavelengths of 
20 nm and 40 nm. Importantly, all of the graphene curves 
show one crossing point with the averaged Si02 curve. 



Raman investigations on monolayer graphene 

In order to observe the influence of the graphene film 
morphology on Raman scattering, we compared Raman 
spectra (laser wave length 532.1 nm), taken in identi- 
cal experimental conditions, of a monolayer flake which 
is partly suspended and one that follows the substrate 
(Fig. 11). The surface structure of these samples was 
determined using STM measurements (Insets Fig. [TT^ i) . 
The Raman spectra were recorded before and after the 
STM measurements giving the same result. A differ- 
ence between both spectra regarding the positions of the 
2D and G peaks is clearly visible. The wave number 
of the 2D peak is 2686.7 cm -1 for the partly suspended 
graphene flake and 2670.5 cm -1 for the flake following 
the substrate. Hence, it results that the change in mor- 
phology causes a shift of the 2D peak by 16.2 cm -1 (Fig. 
[TT^i). A similar effect was observed at the G peak which 
is shifted by 5.5 cm -1 to lower values for the flake fol- 
lowing the substrate (Fig.[TT]3). The D peak (about 1350 
cm -1 ) is barely visible for both samples indicating a good 
quality of our graphene. Raman spectra of several of our 
graphene samples (about 15) prepared in nominally iden- 
tical conditions show more frequently a non shifted 2D 
line with an average value of 2692.4 ±5.3 cm -1 . Hence, 
we conclude that our samples mostly exhibit intrinsic rip- 
pling. If we compare with wave numbers of the 2D line 




FIG. 6: (Color online) AFM images of Si0 2 with tips of different radii: Si tip NSC15 with 10 nm (a) and W/Si tip DP15/HiRes- 
W/A1BS with 1 nm (b), respectively. Both images show the same color scale. Insets: SEM images of similar tips, with kind 
permission of MikroMasch (www.spmtips.com). 
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FIG. 7: (Color online) 2D autocorrelation images of graphene (top row) and S1O2 (bottom row). Left: autocorrelation 
calculated from original images. Middle and right: high-frequency and low- frequency part, obtained by high pass or low pass 
filtering, respectively, using a 20 nm first-order Butterworth filter. Note that the color scale is identical within each column. 
Drift effects have been removed from the graphene images by using the atomic resolution. 



within the literature using only the data taken with a ity of samples investigated so far exhibit intrinsic rippling 
similar laser wave length (514 - 532 nm) jT9l [33j [34], we as displayed in the right inset of Figure 7a. 
find 2680 - 2700 cm -1 . Hence, we believe that the major- 



7 




FIG. 8: (Color online) Different graphene areas (top row), separated by several /mi, with autocorrelation images (middle) and 
high-pass filtered autocorrelation images (bottom row) . The preferential orientation of long-range corrugation varies arbitrarily, 
while the short-range corrugation shows similar orientation. Inset: Autocorrelation of the atomic resolution image. 
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FIG. 11: (Color online) (a) Comparison of Raman spectra at 532.1 nm laser wave length for partly suspended and following 
the substrate graphene monolayer flake. The 2D peak position is 2686.7 cm -1 for partly suspended graphene and 2670.5 cm -1 
for the flake following the substrate. The experimental conditions for both samples are identical. Left inset: STM measurement 
of the graphene flake following the substrate (0.4 V, 0.2 nA). Right inset: STM measurement of the partly suspended graphene 
flake (1 V, 0.2 nA). (b) Raman spectra of the D and G band for the same samples. The D line intensity (1350 cm -1 , marked 
by a dashed line) is very weak for both samples. The G line from the flake following the substrate is shifted to lower values by 
5.5 cm -1 with respect to the G peak of the partly suspended graphene. 



